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{J • ABSTRACT 

X-ray observations of galaxy clusters provide emission measure weighted spectra, aris- 
ing from a range of density and temperature fluctuations in the intra-cluster medium 
(ICM). This is fitted to a single temperature plasma emission model to provide an 
O ', estimate of the gas density and temperature, which are sensitive to the gas inho- 

mogeneities. Therefore, X-ray observations yield a potentially biased estimate of the 
c/2 ' thermal gas pressure, Px- At the same time Sunyaev-Zeldovich (SZ) observations di- 

rectly measure the integrated gas pressure, Psz- If the X-ray pressure profiles are 
strongly biased with respect to the SZ, then one has the possibility to probe the gas 
fNJ \ inhomogeneities (their amplitude and physical nature), even at scales unresolved by 

^ . the current generation of telescopes. At the same time, a weak bias has implications for 

the interchangeable use of mass proxies like Ysz and Yx as cosmological probes. In this 
lO | paper we investigate the dependence of the bias, defined as bp (r) = P x (r) /Psz i r ) — 1 , 

on the characteristics of fluctuations in the ICM taking into account the correlation 
between temperature and density fluctuations. We made a simple prediction of the 
irreducible bias in idealised X-ray vs SZ observations using multi-temperature plasma 
emission model. We also provide a simple fitting form to estimate the bias given the 
CN1 j distribution of fluctuations. In real observations there can be additional complica- 

tions arising from instrumental background, insufficient photon statistics, asphericity, 
method of deprojection, etc. Analysing a sample of 16 clusters extracted from hydro- 
dynamical simulations, we find that the median value of bias is within ±3% within 
■ R500, it decreases to —5% at i? 50 o < r < 1.5 i?5oo and then rises back to ~ 0% at 

r > 2 i?5oo- The scatter of bp(r) between individual relaxed clusters is small - at the 
level of < 0.03 within -R500, but turns significantly larger (0.25) and highly skewed 
(pp(r) 3> 0) at r > 1.5 i?soo- For any relaxed cluster we find |&p(r)| < 15% within 
i?50o, across different implementations of input physics in the simulations. Unrelaxed 
clusters exhibit a larger scatter in bp(r) (both from radius to radius and from cluster 
to cluster). 

Key words: galaxies: clusters: intra-cluster medium — X-rays: galaxies: clusters — 
cosmic background radiation 



43 
Oh 



1 INTRODUCTION 

Galaxy clusters offer an interesting possibility of prob- 
ing physics relevant on both cosm ological as well as 
* satejOmpa-garching.mpg.dc galactic scales (see recent review by iKravtsov fc Borganil 
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12012ft . On the one hand, th ey are useful for c onstrain- 
ing cosmological parameters ([Haiman et al.l 200 if), models 
of dark energy rtMaiumdar fc Mohri l200i ; IVikhlinm" et al.l 
120091 ; [b urcnm fc Vikhlininl l2012f). and possible m odifica- 
tions to the theory of gravity ( Rapetti et al.|[2009h . At the 
same time, they offer insights into feedback processes involv- 
ing the interplay between the intra-cluster medium (ICM) 
and constituent galaxies, and growth of supermassive black 
holes occurring at their centres (e.g.. [Churazov et al.l [20051: 
Croton et all |2006| ; iRaffertv et ail 120061 : ISomerville etail 



2008). 



Galaxy clusters emit strongly in X-rays from the 
shock heating and compres sion of the infalling matter 
(|Sunvaev fc Zeldovichl Il972al ). In the last two decades, a 
large number clusters have been imaged at high resolu- 
tion through X-ray instruments such as ROSAT, XMM- 
Newton and Chandra. The X-ray satellite, Suzaku has 
also been useful at im aging the outskirts of clusters (e.g., 
ISimionescu et alll201ll ), due to its lower instrumental back- 
ground. On the cosmic microwave background (CMB) sky, 
clusters are visible due to the inverse Compton scattering 
of the CMB p hotons by the hot ICM, kno wn commonly as 
the SZ effect (|Sunvaev fc Zeldovichl 1 1972bl). In recen t years 
ongoing SZ surveys li ke SPT (iReichardt et al.ll2012f). ACT 
dMarriage et al.ll201lT ). Planck JPlanck Collaboration et all 
l2011al ) and BOLOCAM jSavers et al.ll2012l ) have been imag- 
ing clusters on arc-minute scales using multi-frequency data. 
These surveys should detect > 1000 clusters up to high red- 
shifts making them valuable as cosmological probes. How- 
ever, in order to place competitive constraints from these ob- 
servations, the cluster mas ses would have to be determined 
with an accuracy of < 5% (|Allen et alj|201ll ). 

Joint observations of clusters in SZ, X-ray and op- 
tical bands provide multiple mass proxies such as inte- 
grated SZ flux, Y S z jMotl et al.l 120051 : lNagalll2006l ); Y x = 
M sli sTx (|Kravtsov et al. 20061 ) from X-ray observations; 



and measurements of velocity dispersions in the optical 
|Evrard et all 120081 ). Forthcoming missions like ASTRO- 
H would probe the gas m o tions via X-ray spectro scopy 
ijlnogamov fc Sunvae y| |2003l : IZhuravleva' et al. I l2012ah and 
help provide unbiased hydro static mass estimates (e.g., 
lLau. Kravtsov. fc NagaH 120091 '). Such multi-wavelength ob- 
servations would be useful to check for the consistency of the 
derived masses and help in identifying the sources of scat- 
ter in the mass-observable scaling relations (Battagli a et al.l 
l2012al : lNoh fc CohrJl2012l ) . Accurate knowledge of this scat- 
ter is important for precise determinati on of the cosmo- 
logical parameters u sing cluster surveys (|Lima fc Hull2005l : 
ICIadders et 111120071 ). 

In this context, some of the recent efforts have focused 
on following up many of the SZ detected clusters in X-rays 
in order to probe the relationship between Yx and Ys z 
l|Andersson et aT | [20lTI : iPlanck Collaboration, et al.1l2011bh , 
both of which measure the total thermal energy in the clus- 
ter. These works have reported the ratio Ysz/Yx at values 
0.82 ± 0.07 and 0.95 ± 0.03 respectively, for me asurements 
within a cluster radius of -RsocQ More recently, iRozo et al.l 



(|2012[ ) independently reported this ratio to be 0.82 ± 0.024 
using cluster data from Planck and Chandra observations. 
More data on joint observations of clusters would reduce 
the statistical and systematic uncertainties and determine 
the bias between Yx and Ysz with better accuracy. 

The deprojection analysis of X -ray surface bright- 
ness, first implemented for clusters bv lFabian et al.l i|l98lf ). 
is be ing routinely performed (in v a rious versions, for 
e.g., Kriss, Cioffi, & Canizaresl Il983l : iMcLaughlirl Il999l : 
IVikhlinin et alj|2006l ) to yield the three-dimensional radial 
profiles of both temperature and density. Henceforth we 
shall use the term 'deproj ection' to imply the particular 
method as implemented inlChurazov et alj (|2003l ) (see also 
iRussell. Sanders, fc Fabianll2008l ). Briefly, assuming spheri- 
cal symmetry, a set of measured values of the surface bright- 
ness in concentric rings is converted into a set of emissivities 
in spherical shells, using the inverse of a simple geometrical 
projection matrix. The procedure is repeated for many en- 
ergies, yielding for each shell an energy dependent volume 
emissivity - i.e. the spectrum, which is then fitted in XSPEC 
with the standard model(s). 

Thanks to rapid technological progress in SZ ob- 
servations it has been possible to apply similar meth- 
ods to well-resolved SZ ima g es obt ained wit h the Planck 
JPlanck Collaboration, et al l |2012| ). SPT (IPlagge et al 
20101) . ACT dReese et al.1 |2012|). CAR MA (jpiagge et a" 



20121). MUSTANG (|Korngut et all 1201 if ) and BOLOCAM 



( Savers et aill2012l ). From these measurements one may then 
obtain the volume integrated thermal pressure J PdV = 
f n e TdV, in each concentric shell both in X-rays and SZ, 
following the deprojection method as outlined in the above 
paragraph. 

X-ray observations reveal the presence of substructures, 
even in relaxed clusters, with the fluctuatio ns in density 
and t emperature having amplit udes of ~ 10% (Zhang et all 
l2009f ) at r < 0.5 7? 5 oo (see also IChurazov et al.M2012f ). Zhu- 
ravleva et al. (2012b), henceforth Z12, investigated prop- 
erties of the ICM fluctuations in simulated clusters and 
found their distribution to be near-lognormal along with a 
high density tail contributing ~ 1% by volume (see also, 
iRoncarelli et al.ll2006l : IVazza et al.ll201lf ). The width of the 
density fluctuations is small (0.1 dex) at the centre of relaxed 
clusters but grows by a factor of 6 at 2 -R500; with unrelaxed 
clusters found to have a broader distribution than relaxed 
clusters. Z12 propose a simple method for identifying the 
high density clumps using a sample of simulated clusters 
and show that the median values of the ICM pressure re- 
main robust to the extraction of such clumps. 

If we assume that the X-ray emissivity is independent 
of temperatur(J3, and the width of the lognormal distribu- 
tion of gas density fluctuations, a ne is small (cr ne <C 1), the 
bias, bp, between the pressures measured from X-ray (Px) 
and SZ (Psz) observations, as a function of a ne , in the fol- 
lowing three different hypotheses (isothermal, adiabatic and 
isobaric gas) may be written as, 



-R500 is the radius, r, of the cluster within which the enclosed 
mass, M(< r), equals, M500 = 500p cr it (z) 3 ; where p C rit(z) is 
the critical density of the universe at the observed redshift, z. 



2 The X-ray emissivity in the energy band E1-E2 is ex oc 
n 2 A(El, E2,T). For T > 3 kcV the soft band (0.5 - 2 keV) emis- 
sivity is almost independent of the temperature. 
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Figure 1. Distribution of density, temperature and pressure in simulated clusters fseel lZl3 . for details). The solid lines show the actual 
distribution, while the dashed lines indicate the lognormal distribution having the same full width at half maxima. The distributions 
shown in blue/red/black colours are measured for the (highly relaxed) cluster CL7 in the CSF simulation (see section[3} in concentric 
shells at the radii, t/R^qq: 0.9—1.0/1.1-1.2/1.6-1.8. Note the presence of an extended high density tail in the distribution of density 
fluctuations. 
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for the adiabatic index, 7 = 5/3. It is interesting to note here 
that the bias takes positive values for isobaric fluctuations 
but is negative for adiabatic fluctuations. 

Thus, for such fluctuations one expects the bias 
to be simply related to the gas clumping factor, 
C = (jig) / (n e ) 2 = exp(a^ e ) ~ 1 + of le . Note also that the 
magnitud e of this bias is small as long as the fluctuations are 
not large. iMathiesen et al.l (1 19991 ) showed using simulations 
that the average bias in gas mass withi n -R500 due to clumpi - 
ness in the ICM is about +16%, while INagai et alJ (12007a ) 
found the bias to be lower at < 6%. INagai fc Laul (|2011 ) 



showed that clumping may be even larger (C ~ 2 — 6) at the 
clust er outskirts (r > 1 — 1.5 R500), (see also IVazza et al.l 

X-ray and SZ observations provide complimentary in- 
formation for studying dumpiness and fluctuations in the 
ICM. However, in reality the actual bias between Px and 
Psz may well be different from the simple estimate pre- 
sented in equation [T] due to the effects of X-ray instrumen- 
tal response and, more importantly, the bias resulting from 
fitting a single temperature model to a spe ctrum produced 
by gas with a range of tem peratures (e.g. iMazzotta et al.l 
12004 iKawahara et alJl2007ft . For example, in a uniform gas 
having only temperature fluctuations, the measured temper- 
ature is weighted down by A(£T, E2, T) due to the relative 



importance of line spectra at lower temperatures, producing 
a negatively biased Px with respect to Psz- In addition, 
there might also be a bias arising from the application of a 
simple ft model (|Rasia et al.| [2006) or from the assumption 
of spherical s ymmetry for cluster s having more complicated 
morphology ( Ameglio et al.|[2~007r ) , however we shall not dis- 
cuss these issues here. 

It is therefore necessary to make a realistic estimate 
of the irreducible bias between X-ray and SZ measurements 
of pressure profiles. By irreducible, we mean the bias that 
would persist even for the most favourable or idealised ob- 
servations. In this work, we estimate this bias in two ways: 
(i) From an idealised hot gas described by a lognormal dis- 
tribution of density and temperature values; (ii) and then 
from hydrodynamical simulations (with different implemen- 
tations of input physics) of clusters in a broad mass range. 
Mock X-ray spectra generated from the given distributions 
of gas density and temperature are used to estimate this 
bias. Our computation of this bias takes into account both 
the density and temperature fluctuations along with their 
correlation, as well as the effects of the X-ray instrumen- 
tal response function (for the Chandra telescope) involved 
in the spectral fitting. We wish to emphasise that in real 
observations, the bias is likely to depend on the method of 
deprojection, X-ray instrumental background, ability to re- 
move dense clumps from images (which may be limited by 
photon statistics), etc. The actual bias inherent in a given 
X-ray/SZ analysis method should be estimated via detailed 
tests using the analysis method outlined in section [5] on the 
mock X-ray observations. 

The remainder of the paper is organised as follows: In 
section[21 we describe the method used to compute the pres- 
sure bias from X-ray and SZ observations and apply it to 
an idealised gas described by a bivariate lognormal distribu- 
tion of density and temperature fluctuations. In section[3]we 
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briefly describe the simulations that were used to estimate 
the bias. In section 0] we discuss the properties of density 
and temperature fluctuations in the ICM of simulated clus- 
ters and define two procedures used to remove the contri- 
bution of high density clumpy regions. Next we present the 
radial profiles of bias in X-ray/SZ pressure measurements 
along with the biases in the observed densities and tem- 
peratures from X-ray observations. In section \S\ we com- 
pare our results with current and future observations and 
also with some recent results presented in other theoretical 
works. Lastly we highlight our main results and conclude. 



2 THE BIAS FROM AN IDEALISED HOT GAS 
WITH A LOGNORMAL DISTRIBUTION 

Hydrodynamical simulations of clusters show that the den- 
sity, n e , and temperat ure, T, of the ICM follow a nearly 
lognormal distribution ()Kawahara et al.ll2007l : l lZ12h with an 
extended tail of high density fluctuations (HDF), see Fig.[T] 
We model an idealised hot gas (T ~ few keV) as a weighted 
sum of 30 x 30 components (in densities and temperature) 
with the weights (here, the emission measures) correctly nor- 
malised to follow a bivariate lognormal distribution in n e 
and T, defined by the following probability density function, 



f(n e ,T) 



■ exp 



\xc~ x x' 1 



(2) 



2im e T\C 

where, X= \hx(n e ) - ix ne , ln(T) - 
and C is the covariance matrix, 

cr„ e 0Wctt£ 

Such a distribution is characterised by the means of the 
lognormal distribution, \x, ne = (ln(n e /l x 10" 
M t = <ln(T/l keV)) (or T mcdian 



i )) and 
exp(pr) keV), the vari- 
ances, a^ e and a^, along with the coefficient of correlation, 
£, between ln(n e ) and ln(T). Of these, one may always set 
the density parameter, n„ e to an arbitrary value since we are 
only interested in computing the ratio Px/Psz', henceforth 
we shall not mention the value of this parameter. 



2.1 Method: Mock X-ray spectrum 

We use the MEKAL code in XSPEC jMewe et all 1 19851 . 
1 19861 ; lKaastralfl99^ ; iLiedahl et ai1ll995h to generate the X- 
ray model spectra from a hot diffuse gas, including the line 
emissions from the astrophysically abundant elements, for 
each of the components defined by the values of n e and T. 
The choice of the MEKAL model was motivated by its abil- 
ity to calculate anew the spectrum at any given tempera- 
ture, rather than interpolate from tabulated values as in the 
APEC modeE 

These spectra are then added together to 
produce a composite spectrum of X-ray emission from an 
ICM having a distribution of temperatures and densities. 
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Figure 2. Single temperature fit to a composite spectrum, for 
the Chandra response. Upper panel: The red curve shows the 
composite spectrum generated from a bivariate lognormal distri- 
bution of gas components with ot = 0.45, cr ne = 0.55 and cor- 
relation £ = —0.55, while the blue curve is the best-fitting single 
temperature MEKAL spectrum. Lower panel: Difference between 
the data and model. For this spectrum the best-fitting values of 
density and temperature are 1.38xl0 — 3 cm -3 and 3.91 keV, giv- 
ing biases of +16.8% and -23.3% w.r.t. n rcf = 1.16 X 10 -3 cm" 3 
and T rc f = 3.91 keV respectively (see equation l3l. Using equa- 
tion [4] the bias in pressures turns out to be -10.4%. A mctal- 
licity value of 0.35 was as sumed relative to the solar abundance 
i Anders fc Grevesselll989h . 



This composite spectrum is next fitted by a single compo- 
nent MEKAL model (in the energy range 0.5 - 8.0 keV) to 
obtain the best-fitting values n e ,flt and Tfn using the Chan- 
dra instrumental response file. The X-ray electron pressure, 
Px = n ei fi t Tfi t , for the ICM is then computed from these 
best-fitting values of density and temperaturfl 

The bias in the gas density is defined with respect to 
the volume-weighted density and the bias in temperature is 
defined with respect to the mass weighted temperature, 



Jn e dV 



T rf ,r = 



J Tn e dV 



JdV ' r " Jn e dV ' 

With these definitions the SZ pressure is, Psz = Pn 
n re fT re f and the bias in pressure may be written as, 

, Px n e at Tm , 
bp = — 1 = — ■ — — 1. 

Psz n e , e f T re f 



(3) 



(4) 



The top panel of Fig. [2] shows the mock X-ray spectrum 
generated from a hot diffuse gas that has a bivariate log- 
normal distribution (red line) with a median temperature, 
Tmedian = 4 keV and median density, n e , mc dian = 1 x 10~ 3 
cm -3 . The blue line shows the best-fitting spectrum spec- 
ified by a single value of temperature and density. Using 
equation|4]we compute the pressure bias, bp — —10.4%. 



3 We find that the APEC model jSmith et alj|200ll) produces 
small spurious jumps in the bias caused by the interpolation of 
the tabulated spectra. These jumps are absent in the MEKAL 
model if the emission spectra are calculated anew for each model. 



4 Note that when fitting a MEKAL model to the composite spec- 
trum we allow only the temperature and normalisation to be free 
parameters, while the redshift, metal abundances, and the Hy- 
drogen column density are kept fixed. 
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Figure 3. Contour map predicting the bias, bp = Px/Psz — li for an idealised hot gas described by a bivariate lognormal distribution 
of density and temperature fluctuations for various values of median temperatures, T me< j; an (in columns) and assumed mctallicities (in 
rows). The bias values are indicated on the colourbar (in percentage) as a function of two parameters a ne and £, using the best fitting 
relation from Fig. [5] cry = 0.73<r r — 0.02. The green contours indicate the zero bias, while adjacent contours are in increments of 2.5%. 
At higher values of T median , the dependence of bias on £ and cr„ e is very similar to that indicated by T median = 3 keV. The bias has 
been computed by fitting a single temperature spectrum in the energy band 0.5-8.0 keV and using the Chandra response function. 



2.2 Results: Bias from a lognormal distribution of 
gas density and temperature fluctuations 

As a preliminary check we verified that our 
scheme reproduces the simple estimate of the ratio 
Px I Psz = n e ,nt /n e ,tet = VC when the temperature fluc- 
tuations are set to zero, while retaining only the density 
perturbations (see also equation [T]). This is done by gen- 
erating a mock spectrum from a lognormal distribution 
of density fluctuations, convolving it with the instrument 
response (Chandra response in our case) and then fitting it 
to a single temperature model, as explained before in section 
I2.ll to compute the bias using equation [4] We then go on to 
systematically explore the dependence of bias in an idealised 
hot gas described by a bivariate lognormal distribution 
of density and temperature fluctuations (see equation [2]). 
For this we compute the bias (as described in section 
I2.ll) as a function of the parameters {ut, One, £, r mcdian } 
for values in the range - 0.01 ^ ^ 0.6, and 

—0.6 ^ £ ^ 0.6 at various temperatures, T mc di an 0.5 keV, 



for an assumed value of metallicity. We find that the bias is 
most sensitive to the value of £, the coefficient of correlation 
between density and temperature fluctuations (see equation 
[2j. Keeping the other parameters fixed, the bias always 
increases almost linearly with £ for |£| < 0.55. The rate of 
increase however depends on the widths of the distribution 
of temperature and density fluctuations, i.e. ot and a ne . 
The simple linear dependence of the bias on £ implies 
the possibility of probing the nature of perturbations in 
the ICM: £ < would indicate the relative importance of 
isobaric or entropj|3 perturbations, while £ > would imply 
that adiabatic or pressure perturbations (caused by sound 
waves or weak shocks) are dominant. 

We now try to describe the dependence of the bias in 
pressure in terms of the parameters of the bivariate log- 
normal distribution using a fitting form. For values of the 

5 Here entropy is defined as S = fc^T '/n^ 3 , where n e and T are 
electron density and temperature of the ICM respectively. 
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Table 1. Values of the best-fit parameters of the fitting form in 
equation \5\ to predict the bias, bp = Px/Psz ~ 1> m X-ray vs SZ 
pressure, valid at temperatures, T mo( ji aI1 ^ 3 keV. 



Z/Zq 


Pi 


f>2 


P3 


0.20 


0.69 


0.75 


i.48 


0.35 


0.75 


0.83 


1.31 


0.50 


0.80 


0.99 


1.08 



median gas temperature, T mc di an ^ 3 keV, the bias may be 
parametrised using three free parameters as follows, 

bp = 7j a ne —p\OT + P2 ln(l + P3&T)cr n e£- (5) 

The best-fitting values of the parameters are indicated in Ta- 
ble[T]for a range of assumed metallicites. This fitting formula 
is accurate (in absolute terms) to within 0.01 for <7 ne / T ^ 0.3 
with |£| ^ 0.6 and within 0.05 for <J ne / T ^ 0.6 with the same 
range for £. We verified that these conditions hold for the 
ICM within Rsoo, in simulated clusters (see Section [4T3] and 
the discussion on Fig. [7)l . For median temperatures greater 
than 3 keV, the X-ray emissivity in the 0.5-8 keV band is 
only a weak function of temperature; at lower temperatures 
however, there is a significant dependence due to relative 
importance of discrete spectral line emission, causing the 
measured bias to depart significantly from this parametri- 
sation. We find that a decrease in the median temperature 
below ~ 1 keV causes the bias to change steeply in a non- 
trivial manner, which is difficult to parametrise through a 
simple form. 

In general the bias in pressure for an assumed value of 
metal abundances depends on 4 parameters £ and 

Tmedian- However, we may simplify this by eliminating a pa- 
rameter using an approximate relation, <tt = 0.73<tt — 0.02, 
that is seen in hydrodynamical simulations of clusters (see 
Fig. [5|. With this simplification, in Fig. [3] we plot the bias 
contours for various values of median temperatures as a 
function of £ and a ne assuming two values of gas metal- 
licities, Z/Zq — 0.2 and 0.5; the former being the typ- 
ical value beyond 0.2 Riso, while the latter is indicative 
of the ICM metallicities at the centre of the cluster (e.g. 
iLeccardi fc Molendilliofih . We see that for ICM tempera- 
tures, Tmcdian > 3 keV, the bias is small (from —10% to 
+5%) for most of the parameter space in the (£, a n e) plane. 
Only for higher values of both the fluctuations and the neg- 
ative correlation £, the bias goes beyond —15%. For lower 
temperatures (1 keV < T mc di an < 3 keV), this bias can be be- 
tween —25% to +15%, depending on the metal abundance; 
while for even lower temperatures the bias turns positive, 
mostly lying at the — 10% level. Note that especially for 
low temperatures (T med ia n < 3 keV) the value of bias is more 
sensitive to the assumed metal abundance; the effect of in- 
creasing metallicity is to reduce the bias when positive, or 
to move it towards more negative values otherwise. 



3 SIMULATIONS AND SAMPLE OF GALAXY 
CLUSTERS 

We use a sample of 16 simulated clusters at z — taken 
from iNagai et all (|2007al lbl) . The simulations were done us- 
ing the Adaptive Refinement Tree (ART) N-body+gas- 



Table 2. Propert ies of simulated clusters in the sample of 
INagai et al, (2007a b) at z = 0. Two runs were performed with 
different physics, CSF run with cooling+star formation and NR 
run without radiative cooling and star formation. The last col- 
umn indicates if the cluster is visually identified as relaxed (R) 
or unrelaxed (U). 

Cluster ID R 500 M'g* Relaxed (R) or 
(fe- 1 Mpc) (10 14 fe- 1 M Q ) Unrelaxed (U) 



CSF / NR CSF / NR CSF / NR 



CL101 


1.16 


/ 


1.14 


9.08 


/ 


8.62 


U / U 


CL102 


0.98 


/ 
/ 


0.95 


5.45 


/ 


4.63 


u / u 


CL103 


0.99 


/ 


0.99 


5.71 


/ 


5.71 


u / u 


CL104 


0.97 


/ 


0.97 


5.39 


/ 


5.31 


R / R 


CL105 


0.94 


/ 


0.92 


4.86 


/ 


4.50 


u / u 


CL106 


0.84 


/ 


0.84 


3.47 


/ 


3.40 


u / u 


CL107 


0.76 


/ 


0.78 


2.57 


/ 


2.74 


u / u 


CL3 


0.71 


/ 


0.70 


2.09 


/ 


1.98 


R / R 


CL5 


0.61 


/ 


0.61 


1.31 


/ 


1.34 


R / U 


CL6 


0.66 


/ 


0.61 


1.68 


/ 


1.32 


U / R 


CL7 


0.62 


/ 


0.60 


1.41 


/ 


1.25 


R / R 


CL9 


0.52 


/ 


0.51 


0.823 


/ 


0.775 


u / u 


CLIO 


0.49 


/ 


0.47 


0.672 


/ 


0.621 


R / R 


CL11 


0.54 


/ 


0.44 


0.899 


/ 


0.482 


U / R 


CL14 


0.51 


/ 


0.48 


0.769 


/ 


0.652 


R / R 


CL24 


0.39 


/ 


0.39 


0.347 


/ 


0.347 


u / u 



dynamics code (jKravtsov et al ] |l997l . 120021 ) assuming a flat 
ACDM cosmology with the following values for the cosmo- 
logical parameters: Q m = 0.3, Qt = 0.04286, h = 0.7 and 
erg = 0.9. Two types of simulations were used with differ- 
ent physics involved, but with the same initial conditions: 
(i) non-radiative (NR) runs without any radiative cooling 
and (ii) cooling plus star formation (CSF) runs, which in- 
cluded metallicity-dependent radiative cooling, star forma- 
tion, supernova feedback and UV background. The clusters 
in the sample were selected randomly to uniformly sample 
the mass range, 7xl0 13 h" 1 Mq < M 50 o < 2xl0 15 h' 1 M . 
These clusters were selected from low-resolution simulations 
and resimulated at a higher resolution. The division of the 
sample into relaxed and unrelaxe d sub-samples (see Table[2} 
was done in INagai et alj (|2007al ) by visually examining the 
morphology of mock X-ray images in three different projec- 
tions. This was done by producing mock images of the sim- 
ulated clusters using the Chandra response function with 
a 100 ks exposure. From these images the 'relaxed' clusters 
were identified as those with a regular X-ray morphology and 
without any significant deviations from elliptical symmetry. 
On the other hand, clusters were classified as 'unrelaxed' 
when departures from elliptical symmetry, filamentary X- 
ray structures, or large isophotal centroid shifts were seen 
in the produced images. This procedure is similar to that 
used by observers. Z12 made another classification based on 
the widths of the density distribution, er ne , at r = R500. This 
classific ation agrees very well with that done in lNagai et al.l 
(|2007af ). and the difference in CT„ e (i?5oo) between relaxed 
and unrelaxed clusters was at the level of m 2.5a. There 
also exist other schemes of classification, for e.g. using the 
two dimensional m ultipole expansion of the projected grav- 
itational poten tial (iBuote fc Tsail Il995l ) or computing the 
centroid shifts (|Poole et al.ll2006r i7 however we did not com- 
pare our classification with these. 

The ART code uses a non-uniform and dynamic mesh 
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adjusting to the evolving particle distribution, thus allowing 
for a large dynamic range. To simplify our analysis, the code 
output (gas densities, temperatures, velocities, and metal- 
licities) was interpolated onto a uniform grid with grid cells 
of size comparable to the effective resolution of the clus- 
ter simulations. Following Z12, we randomly sample 4 x 10 7 
data points within the sphere of radius of 5 h' 1 Mpc with 
a weight oc 1/r 2 , where r is the distance from the cluster 
centre. This sampling is uniform in azimuthal and polar di- 
rections and provides equal number of points per spherical 
shell of a given thickness. Such a scheme is useful for dealing 
with averaged quantities computed over spherical shells as 
discussed further in the following sections. 



4 THE BIAS FROM SIMULATED GALAXY 
CLUSTERS 

In this section we shall compute the bias between 'measured' 
pressure profiles from X-ray and SZ observations using the 
sample of simulated clusters introduced in the previous sec- 
tion. We wish to emphasise that the term 'measured' (here 
and henceforth) refers to the value of pressure as would be 
inferred from observations in X-ray/SZ through a deprojec- 
tion analysis of real data. Note that this may well be different 
from the actual values of thermal pressure of the ICM. 

Each cluster is divided into 32 concentric shells (spaced 
in equal logarithmic intervals) between the radii 0.1 < 
r/Rmo < 2.3, in which we compute the histogram of den- 
sity and temperature values. This is then used to generate 
fake spectra assuming a metallicity of Z/Zq = 0.35, as 
described before in section 12.11 The fake spectra are then 
analysed using the Chandra response matrix to give the 
best-fitting values of density, n e ,at and temperature, Tat 
for each shell. The 'measured' X-ray pressure is then sim- 
ply Px = ra e ,fitlfif On the other hand the SZ pressure for 
each shell is obtained by computing the mean thermal pres- 
sure in the shell Psz = Yli n e,iTiWi/ w ii where Wi is 
the weight, i.e. the fractional volume occupied by the ICM 
specified by the values n ei i and Ti for density and tempera- 
ture respectively. These are used now to compute the bias, 
bp(r) = Px(r)/Psz — 1 between the 'measured' X-ray and 
SZ pressures, in each shell. 

4.1 Properties of the ICM fluctuations 

Z12 used the same set of simulations, as described in sec- 
tion [21 to study the properties of the inhomogeneities in 
the ICM. They developed a simple and robust method to 
separate the HDF from the nearly hydrostatic bulk compo- 
nent. This was implemented by imposing a cut within each 
shell to exclude all particles having densities, n e , such that 
log(n e ) > log(n e , mc dian) + 3.5<7ne- Henceforth we shall use 

6 Strictly speaking, assuming a fixed metallicity value for differ- 
ent shells is reasonable only within ~ 0.5 -R500 where T mc( ji a n > 
1 keV for most of the simulated clusters, and the bias depends 
only weakly on the assumed metal abundance. For T mc( j; an < 1 
keV, we find that our results on bias are more influenced by metal- 
licity; when higher accuracy is required, especially so for lower 
temperatures, one has to do the direct analysis on a case by case 
basis. 



the term 'HDF' to exclusively denote the particles in the 
simulation lying above this cut, and the term 'bulk' to de- 
scribe the particles lying within this cut. 

Here we review some of the results presented in Z12. 
Fig. U shows the radial profiles indicating the properties of 
ICM fluctuations for the sample of relaxed clusters in the NR 
simulations. In Fig. [4] (a) the solid lines denote the median 
values (over the sample of relaxed clusters in the NR run) 
of the parameters obtained by fitting a bivariate lognormal 
distribution (see equation [2j within each shell. The dashed 
lines are the same quantities computed after the removal of 
the HDF components. The thin lines indicate the maximum 
and minimum outliers of the sample that shows the scatter 
across the sample. Fig. [5] shows the relation between the 
parameters a ne and ctt in each shell; the straight line is the 
best-fitting linear relation, ot = 0.73<r, le — 0.02, which is also 
the assumption used in producing the contours in Fig. [3] 
The data from NR and CSF simulations indicate an intrinsic 
scatter of 0.07 in gt for a given value of a ne in this relation. 
Fig- H] (a) and (b) show that the values of the parameters, 
<T„ e and £, are fairly robust to the exclusion of the HDF 
components from the bulk. 

In both the CSF and NR runs, the unrelaxed clusters 
have not only higher values of temperature and density fluc- 
tuations, but also show significantly larger scatter from clus- 
ter to cluster. Both relaxed and unrelaxed clusters show a 
steady increase in the fluctuation amplitude with radius, 
which steepens in the outskirts (r > Rsoo)- As shown in 
Fig- S] (a), the fluctuations in gas density, a„ e , are large 
beyond -R500, and these fluctuations are reduced if the HDF 
components are removed. Fig. [4] (b) shows that the corre- 
lation between temperature and density fluctuations, £, is 
mostly negative, but increases towards zero at larger radii. 
The increase in £ on removal of high density clumps is in 
line with our expectations, as the clumpy regions contribute 
to an anti-correlation (or £ < 0) between density and tem- 
perature fluctuations. 

4.2 Removing the contribution from the brightest 
clumpy regions in the cluster outskirts 

We expect the 'measured' quantities related to the ICM 
(e.g., gas pressure, density and temperature) to be less bi- 
ased when the HDF components are removed from the sim- 
ulation data. Clearly the removal of HDF according to the 
criterion described in the previous section is only possible 
in simulations. In real observations, the complete removal of 
HDF seems unlikely, due to the presence of the X-ray back- 
ground noise and limited photon statistics, especially rele- 
vant at the cluster outskirts or for the high redshift clusters. 
Here we shall attempt to check this by dealing directly with 
the 2-dimensional projected X-ray images obtained from the 
simulated clusters, only in the limit of infinite photon statis- 
tics (i.e., without any Poisson noise) and in the absence of 
instrumental background. However, we shall argue in section 
14. 3 1 that even a partial removal of the clumps, especially the 
brightest ones (as permitted by the photon statistics) , helps 
to reduce the bias to a large extent. 

For each of the clusters, we divide the X-ray images 
into six radially equal annular regions starting from r = 
-R500 as the innermost ring and up to r = 2.3 -Rsoo- We 
choose 7?soo as the inner radius for the purpose of remov- 
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Figure 5. The coloured dots indicate the values of fluctuations 
in density, cr ne , and temperature, cry, in concentric shells from 
the combined sample of clusters in both the NR and CSF simu- 
lations. The colours (mapped to colourbar on right) indicate the 
radial distance of the shell from the cluster centre. The black line 
indicates the best fitting linear relation, cry = 0.73cr ne —0.02. For 
this relation the scatter in cry for a given value of cr ne is 0.07. 



ing clumps since we find that the bias outliers become 
large only beyond this radius, suggesting that the effect 
of clumps is dominant in the bias p rofiles only in the out- 
skirts (r > R500) of clusters (see also [Mathiesen et al.| [l99£ 



iNagai fc Laul201ll ; iBattaglia et~aill2012bl ; IVazza et aljboiS , 
and Fig. 0. Using the pixels lying within each annular 
region we then create histograms of the X-ray brightness 
(in the soft band from 0.5-2.0 keV using the Chandra re- 
sponse) and mask all pixels with brightness, B, such that 



log(B) > log(B mcd ian) + /cutcs, to zero values, where a B is 
the standard deviation of the logarithm of pixel brightness 
in each annulus. Using three different values of / cut - 2, 2.5 
and 3, we then eliminate all the particles from the simu- 
lation data which lie exactly along the line of sight of the 
masked pixels, some of which are responsible for producing 
the very bright regions in the outskirts of the X-ray image. 
We shall henceforth refer to this method of removing the 
high-density contributions from the X-ray images as the 2- 
D cut in contrast with the 3-D cut performed directly on 
the radial shells and as outlined in section [4~T1 and detailed 
in section 4 of Z12. 

We note that the 2-D cuts do not remove all of the 
clumpy contributions from the 3-dimensional radial shells, 
partly because clumps lying along the line of sight through 
the bright central region of the projected image of the clus- 
ter would not show up with high enough contrast in the 
brightness images, and would remain in the filtered (i.e. af- 
ter a 2-D cut) simulation data. For a fair comparison of the 
results from 2-D cut and 3-D cut we attempt to correct this 
by applying a projection in another plane and repeating the 
masking procedure in this new plane. Fig. [6] shows a com- 
parison of projected X-ray images of the simulated cluster 
CL101 under 2-D (using / cut = 2.5) and 3-D cuts. On the 
left is the original image produced without any cuts, while 
the central and the rightmost images are after application 
of 3-D and 2-D cuts respectively. Note that the 3-D cut is 
applied everywhere while the 2-D cut is only applied within 
-R500 < r < 2.3 -R500 (indicated by the annulus which is 
bounded by green circles). We see that the bright clumps 
in this annulus, visible prominently on the leftmost image 
are absent after applying either the 3-D or 2-D cuts. The 
fact that this annulus appears very similar after 2-D and 
3-D cuts indicates that both of these cuts remove nearly the 
same points from the 3-D shells. 

We now try to examine if a similar 2-D cut (or mask- 
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.00e+00 1.00e-09 3.01 e-09 7.04e-09 1.50e-08 3.12e-08 6.30e-08 1.27e-07 2.55e-07 5.08e-07 1.01e-0 

Figure 6. Exclusion of the bright clumps from the projected soft band (0.5—2.0 keV) X-ray image of the simulated cluster CL101 in 
the NR run. The images from left to right are: (i) Original image; (ii) Image showing only the bulk component (i.e. with 3-D cut); 
(iii) Image after masking of the brightest regions (i.e. with 2-D cut specified by the value /cut = 2.5). The 2-D mask is applied in the 
region -R500 < r < 2.3 -R500 whose boundaries are indicated by the green circles, while the 3-D mask is applied everywhere. 



ing) is possible in real observations, where the ability to 
mask the clumps would depend on the signal to noise ratio 
at r > -Rsoo- To be able to identify and mask such regions it 
is necessary that the intrinsic fluctuations of surface bright- 
ness should dominate over the Poisson fluctuations of the 
photon counts. We now estimate the minimum number of 
counts required in order to be able to mask such clumps, 
provided that they are resolved by the instrument. Let N 
be the total number of counts, n rcg be the number of inde- 
pendent regions (or pixels used for masking); then, we have 
n = iV/nrcg as the number of counts per region. To reliably 
remove the region with the surface brightness X times the 
mean value (i.e. Xn counts) one needs (case with no back- 
ground) as a conservative estimate: (X — l)n > 5y/n for 
n >> 1, or (X — l)n > 5 for n « 1. For the case shown in 
Fig. [6] we use 63 by 63 kpc pixels, yielding iV rC g = 9048 in 
the annulus from R500 to 2.3 -Rsoo- Assuming 100,000 counts 
for the whole image, and 7.0% flux coming from this annu- 
lus, the typical number of counts per pixel, N/n las = 0.78. 
Thus, pixels with the surface brightness X > 6.4 are easily 
detectable as more than 5a deviations on top of the mean 
level. To get 100,000 photons in the 0.5-2.0 keV band the 
Chandra ACIS-I instrument with an effective are 43 of ~370 
cm 2 would require an exposure time of ~10 ks for the sim- 
ulated cluster CL101 with a mass of ~ 1.2 x 10 15 M Q and 
~150 ks for CL24 with a mass of ~ 5.0 x 10 13 M , assuming 
a redshift of z = 0.1 for both objects. Further, we see that 
the clumps in the image will still be resolved if the size of 
the pixels is increased by a factor of ~ 2.5. This reduces the 
required number of photons, we now find that pixels with 
surface brightness X > 10.25 can be used for masking, using 



This is the weighted average effective area over the energy range 
0.5—2.0 keV for a hot gas with T ~ few keV. The photon rate used 
to estimate the exposure was also computed for the same energy 
range. 



/cut = 2.5 and 3.0, even if there are 10,000 counts for the 
whole image. 



4.3 Results: Bias between X-ray and SZ pressure 
profiles 

Fig. [TJa) shows the median bias profile of the relaxed clus- 
ters in the NR run. The profiles with and without HDF com- 
ponents (the 3-D cut) are shown by the thick solid blue and 
thick dashed magenta lines, respectively. Thick dot-dashed 
black line shows the median profile with the exclusion of the 
bright clumps beyond R500 in the projected map (the 2-D 
cut with /cut = 2.5). The sample outliers are indicated by 
the thin lines of similar type. Note that this bias is com- 
puted from the composite spectra corresponding to the ac- 
tual distribution of temperature and density values within 
each shell (see section [2~T]| . and not from the idealised log- 
normal distribution. We find the median bias computed from 
the relaxed clusters to be within 0% to +2% at r < 0.5 Rsoo- 
This then decreases to negative values ~ —6% at r = Rsoo- 
The bias profiles remain unaffected by the HDF component 
(due to small fluctuations) within R500, producing almost 
identical profiles with and without HDF. Beyond R500, the 
median bias profiles from 'only bulk' follow the same trend 
as 'bulk+HDF', but are smoother and are limited to within 
—2.5% to +2%. With the 2-D cuts, the median bias profile 
is also smoother and agrees with the bias profile after the 
3-D cut (see the lower panel of Fig. [SJ). 

We find that the unrelaxed clusters display a larger neg- 
ative bias and larger scatter, both from shell to shell and 
across the sample at a given value of r/R 5 oo- The cluster-to- 
cluster scatter increases with radius. At r = -Rsoo, the scatter 
is less than 3% and 13% for relaxed and unrelaxed clusters, 
respectively, for the NR run. At larger radii (r ~ 2 R500), 
the scatter is highly skewed in the positive direction, with 
some shells showing dp>0 (see Fig. [SJ , from the presence of 
dense clumps. After the 3-D cut, the scatter in the pressure 
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Figure 7. Radial profiles showing various biases computed in concentric radial shells of a sample of relaxed clusters in the non-radiative 
simulations. In all the figures, the thick lines indicate the median profile while the thin lines are the most extreme outliers in the sample. 

(a) Pressure bias, bp = Px/Psz — lj for the 'measured' X-ray (Px) and SZ (Psz) pressures. 

(b) The ratio of the 'measured' X-ray/SZ pressures to the median pressure shown in thick red/blue lines for 'bulk+HDF'. 

(c) Bias in the 'measured' density from X-ray observations with respect to n rG f (see equation [3j . 

(d) Bias in the 'measured' temperature with respect to the mass weighted temperature, T re f (see equation[3j. The 'measured' temperature 
is determined simultaneously with the best- fitting n e . 



bias drops to less than 7% within -R500 even in unrelaxed 
clusters. 

Fig. [TJb) shows the 'measured' values of the X-ray and 
SZ pressures w.r.t. the median pressure in concentric shells 
at various radii. Since the value of Psz actually corresponds 
to the average electron pressure, the fact that there is a 
bias even in the SZ profiles implies that the mean and the 
median values of pressure fluctuations are different. This 
is true for any asymmetric distribution such as the near- 
lognormal distribution of pressure fluctuations seen in Fig. [I] 
The plot shows that the 'measured' SZ and X-ray pressures 
are largely similar, the differences between them are smaller 
than differences in the mean and median values of pressure 
fluctuations. 

In Fig. E|c) and (d), we also plot the bias in the den- 
sity and spectroscopic temperature values as would be 'mea- 
sured' from X-ray observations. As seen from equation [T] 



the best-fitting values of n e are biased in the positive direc- 
tion, while the best-fitting Tx is biased negatively. For the 
relaxed clusters in the NR simulation, the median bias in 
gas density increases slowly to +3.5% at -R500 but can reach 
~ +40% at larger radii. Interestingly, the exclusion of bright 
clumps reduces the median bias to +6.5% at r — 2 i?soo 
and makes the bias profile considerably smoother. Beyond 
-R500 the median bias profiles in gas density obtained after 
both 3-D and 2-D cuts are quite similar. The median bias in 
temperature takes similar values up to r ~ 0.8 -R500 as the 
density bias, but in the negative direction, reaching ~ —14% 
at r = 2 i?5oo- We find that the median bias in temperature 
reduces only slightly even after removing the HDF, because 
the removal of the HDF components only affects the den- 
sity distribution, but does not directly limit the temperature 
fluctuations. This is because, there is almost zero correlation 
between density and temperature fluctuations (see Fig. 2] 
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(b)) in the ICM at r > -Rsoo- At the same time, the sample 
outliers producing a negative bias in temperature reduce to 
nearly a third of the original value after the 2-D and 3-D 
cuts (as the very dense clumps are cold). 

In Fig. [8] we compare the bias profiles in pressure, in 
the cluster outskirts (r > -R500), for the various cuts in sur- 
face brightness applied in concentric annuli on the projected 
images of the NR simulations. We also plot for comparison 
the results from the 'HDF+bulk' (No cut) and 'only bulk' 
(3-D cut). We find that application of more stringent cuts 
(specified by lower values of /cut) reduces the scatter in the 
bias profiles. Exclusion of the bright clumps in the clus- 
ter outskirts (r > -R500) reduces the outliers (on the posi- 
tive side) in the bias (bp) profiles, by more than an order 
of magnitude, from +8.2 to +1.8, +1.2, and +0.7 respec- 
tively, at 2.2 i? 5 (M, for / cut = 3.0, 2.5 and 2.0. For / cut = 2.0 
the excluded clumps have surface brightness between ~ 9 
- 13 times that of the average value (across the inner to 
the outer annuli). For / cut = 2.5 this is 13-18; while for 
/cut = 3.0 these numbers are 20-28. It is clear that such ex- 
clusions would remove only the brightest regions associated 
with high density clumps, as shown in Fig. [5] for the unre- 
lated cluster CL101 in the NR run. On the negative side 
(i.e. outliers with bp < 0), the various cuts do not affect the 
outliers as much as for the positive side. We see from the 
bias profiles plotted in the upper panel of Fig. [8] that the 
2-D cut specified by the value / cu t = 2.5 produces similar 
results to the 3-D cut. A more stringent cut (/cut = 2.0) 
may further eliminate the intrinsic brightness fluctuations 
associated with inhomogeneities in the diffuse ICM, while a 
more relaxed cut (/ cu t = 3.0) does not remove all the visibly 
bright clumps. In the lower panel of Fig. [8] we indicate the 
median values of the bias for the 16 clusters; these are again 
quite robust to the application of the various 2-D and 3-D 
cuts. 



Finally, we compare the bias in the median pressure 
profiles, bp = Px/Psz — 1, for both relaxed and unrelaxed 
clusters in simulations with very different input physics. As 
expected, we find more clumping in the CSF simulations 
when compared to the NR simulations; the bias profiles in 
the CSF runs also display significantly larger scatter from 
shell to shell (over the same radial range) due to the presence 
of some very dense clumps. Fig. [9] shows the median bias 
profiles to be within —4% to +2% for r < 0.5 -R500 and 
within —6% to +2% for 0.5 R500 < r < R500 for relaxed 
and unrelaxed clusters in both NR and CSF simulations. 
At r ~ 1.2 -R500, the bias reaches —13%, arising from the 
presence of infalling clumps, but reduces to less than —5% 
when the HDF component is removed. In both the NR and 
CSF simulations the bias increases beyond r > 2 Rsoo- 

More importantly, we find that the pressure bias is lim- 
ited to within ±15% for any of the relaxed simulated clusters 
in both CSF and NR runs of the simulation inside ifeoo- 
Also, the results on median bias profiles are quite similar 
(see Fig. O across very different input physics used in the 
NR and CSF runs. This suggests that our results are robust 
to possible uncertainties in the implementation of physical 
processes in the hydro codes. 



1.2 



1.4 



r/R 5( 

1.6 



2.2 



-No cut (bulk+HDF) 

- 3a cut 

- 2.5a cut 

- 2a cut 

- only bulk 




Figure 8. Reduction of the bias, bp = Px/Psz — 1, in the 
pressure profiles by the introduction of various cuts in surface 
brightness (applied to the X-ray images), at r ^ Rsoo- Upper 
panel: Maximal deviations (or sample outliers) for both relaxed 
and unrelaxed clusters in the NR simulation. Lower panel: Me- 
dian bias profiles for each of the above cuts. Note that in both 
panels, the left axis shows bp = Px / Psz ~ 1 while the right axis 
displays P x /Psz\ the dotted line, 'ZERO BIAS' indicates b P = 
or Px/Psz = 1. 



5 COMPARISON WITH CURRENT AND 
FUTURE MEASUREMENTS 

High-resolution observations of galaxy clusters jointly in X- 
rays and millimetre/sub-millimetre bands provide indepen- 
dent measures of the thermal pressure profiles through a 
deprojection analysis. In this work we provide an estimate 
of the expected bias in the pressure profiles derived from 
the X-ray and microwave bands taking into account both 
temperature and density fluctuations in the ICM, as well 
as the correlation between them. We also predict biases 
in three-dimensional density and temperature structure de- 
rived from X-ray observations. The bias in density depends 
on the ability (determined by the photon statistics and X- 
ray background) to mask out X-ray emitting gas clumps in 
the outskirts of galaxy clusters. Our estimates from section 
4.2l seem to suggest that X-ray observations would be able to 
identify these clumps (if resolved) in clusters up to moderate 
redshifts (z ~ 0.1). 

Accurate measurements of the bias as a function of clus- 
ter radius will allow us to probe the inhomogeneity and 
the nature of perturbations (pressure vs. entropy) in the 
ICM. Our work is also relevant for the use of galaxy clus- 
ters in cosmological studies. A robust characterisation of 
biases between X-ray and SZ observations is critical for con- 
straining cosmological parameters through measurements 
of the angular diameter distances, cIa, at various red shifts 
i|Bonamente et al.ll2006l ; iKhedekar fc Maiumdarjl2010t ). The 
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Figure 9. The bias, bp(r) = Px/Psz — lj in pressure profiles 
for a sample of 16 clusters (both relaxed and unrelaxed) in the 
NR and CSF simulations. For the figure in lower panel, labelled 
as 'only bulk', the bias is computed after excluding contributions 
from the high density fluctuations (HDF) associated with dense 
clumps. 



prediction of the bias in density profiles is also important 
for the determination of the baryon budg et in clusters (e.g., 
ISimionescu et al.|[201ll : IVazza et al]|2012l ). 

Our results appear to be consistent with the bias 
(—10% to —40%, from the cluster centre to r = 
0.8 -R500) in the pressure profiles seen in the Coma cluster 
(jPlanck Collaboration, et al.ll20ll ). using data from XMM- 
Newton and Planck observations. Coma is an unrelaxed clus- 
ter and the bias outliers seen in the small sample of un- 
relaxed simulated clusters range from —15% to +20% and 
-60% to +20%, for the NR and CSF simulations respec- 
tively, in the same radial range. The current results on the 
ratio Ysz/Yx appear to be mixed and inconsistent with 
each othe r. On the one hand, rec ent observations from the 
CARMA l|Bonamente et al.ll2012h indicate that this ratio is 
consistent with unity, while other mea surements (for e.g., 
iRozo et al.ll2012l ; lAndersson et al.ll201lf ) are quite different 
from the biases derived from our simulations. In our view 
this is likely due to some underlying systematic uncertain- 
ties associated with the interpretatio n of the current ob- 
servations (see iBonamente et alj [20061 . for a list of possible 
systematic err ors). 

Recently, iBattaglia et~aH (|2012bl ) estimated biases in 
the gas mass fraction, / gas , derived from X-ray observations, 
using a set of cosmological cluster simulations. Our results 
on the bias in the derived gas density from X-ray data due to 
dumpiness are consistent with their findings. They estimate 



the contribution of clumping to Af gas (< -R200) to be 10 — 
—20% while our results indicate a median bias of 16% in the 
density at the same radius. 



6 CONCLUSIONS 

The main conclusions of this paper are summarised below: 

• We first investigate the dependence of bias between the 
X-ray and SZ pressures by fitting the X-ray emission from 
an idealised hot gas having a lognormal distribution of den- 
sity and temperature values. Our results demonstrate that 
the bias is expected to be small (within ±10%), as long 
as the magnitude of density and temperature fluctuations 
is sufficiently small (err, u ne < 0.45) and their correlation 
does not take large negative values (£ > —0.3), for median 
temperatures, T mc di an ^ 3 keV. However, the bias may be 
somewhat larger (~ 15%) for clusters with lower tempera- 
tures (T me dian < 1 keV) . We provide a fitting form to predict 
the bias as a function of the properties of fluctuations of the 
ICM for T me dian ^ 3 keV and plot bias contours in the (£, 
a) plane that may be used to probe the properties of fluc- 
tuations from a given observed ratio of Px/Psz- At lower 
temperatures (Tmedian < 3 keV) , the bias is also found to be 
a function of the metal abundance in the gas. 

• We then use a sample of 16 simulated clusters to study 
properties of the ICM fluctuations. The clusters in the sim- 
ulations show an increase in the amplitude of both density 
and temperature fluctuations with radius. The fluctuations 
of temperature and density are mostly negatively correlated 
within r < J?soo- This correlation is close to zero at higher 
radii. We compute Px and Psz profiles using the ICM from 
simulated clusters as would be measured from X-ray and SZ 
observations. Both Psz and Px are biased with respect to 
the median pressures values. For Psz this bias is purely due 
to the asymmetric (near lognormal) distribution of pressure 
fluctuations; for Px , in addition to the bias from asymmetry, 
there is also a bias arising from fitting a single temperature 
model to the multi-component emission spectrum. 

• We show that the median bias, bp(r) = 
Px(r)/Psz(r) — 1, between the pressure estimated 
from X-ray and SZ observations is small and lies within 
-6% to +2%, up to -R500, even for unrelaxed clusters. The 
scatter in the pressure bias is significantly smaller for 
relaxed clusters (< 0.03) than unrelaxed clusters (< 0.13) 
at r < -R500; however, it becomes large and positively 
skewed at r ~ 2 i?5oo- For non-relaxed clusters, there is a 
noticeably higher scatter in the bias values from shell to 
shell. The outliers responsible for the scatter in bias can 
be reduced by almost an order of magnitude just by using 
a simple method of excluding the brightest clumps (as 
permitted by reasonably good photon statistics) from the 
cluster outskirts (r > -R500), for both relaxed and unrelaxed 
clusters. 

• Finally, our results on the median values of pressure bias 
between X-ray and SZ observations are not significantly dif- 
ferent in the two types of simulations with different input 
physics: (i) Non-radiative (NR) runs and (ii) Cooling + star 
formation (CSF) runs. Therefore these results, (i.e. a small 
bias within -R500) may be considered to be robust with re- 
spect to possible uncertainties in the physical implementa- 
tions in the hydro codes. 
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